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Abstract. For configurational space of arbitrary dimension a strict form of the 
uncertainty principle has been obtained, which takes into account the dependence 
' of inequality limit on the effective number of pure states present in given statistical 

iy) . mixture. It is shown that in a state with minimal uncertainty the density operators 

eigenfunctions coincide with the stationary wavefunctions of a multidimensional 
harmonic oscillator. The mixed state obtained has a permutational symmetry which 
is typical for a system of identical bosons. 
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1. Introduction 

As shown in the papers [1-3], the generalization of conventional uncertainty relation 
of the position-momentum type (see, e. g. Ref. 4) to the case of mixed states, 
determined by a density matrix p(x,x'), leads to a radically new result. It sounds 
as follows — the smallest allowable room Av p , occupied by system in the phase space, 
grows proportionally to the effective number N e g of pure states by whose statistical 
mixture the given density matrix is representable, that is 

min(A?7p) oc N eff . (1) 

One of the possible variants of relation (1), which can be used in a real or configurational 
space X = (xi, . . . , x s ) with an arbitrary number of dimensions s, has been obtained in 
[5] and is of the form 
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The next definitions are introduced here: Ax, Aq the root mean square widths 
(additionally averaged over all s degrees of freedom) correspondingly of coordinate and 
wave-vector distributions in the state p(x, x') 



(Ax) 2 = - I d s IXV(X,X), (3) 



s 



(Ag) 2 = -i J al s X J d s X'5(X-X') V^p(X,X'), 



2 

Xi 

i=l 



(4) 



N eff = (^J d s X J d s X' |p(X,X')| 2 ) . (5) 
Besides, it is assumed that the normalisation condition 

cfXp(X,X) = l, (6) 



is met and the mean values of position and momentum are equal to zero (without loss 
of generality the latter can be done by choosing an appropriate reference frame). For 
convenience the right-hand side of inequality (2) is represented as two factors. One 
of which, (l/2) s , is the minimum of usual uncertainty relation for pure states and, as 
known, is attained for Gaussian wave-packets [4]. The second factor, C(s), can be named 
the "packing coefficient" and expresses the specificity of N e jj definition and, properly, 
the influence of space multidimensionality [5] 

The particular interest for the relation (2) is caused by its analogy with one of the basic 
statements of statistical mechanics [6, 7] concerning the partitioning of the phase space 
into cells with each of them corresponds to one quantum state. Indeed, the quantity 
Av p = (Ax, Aq) s is a peculiar kind of measure of a system phase volume, whereas N e g 
is a characteristic number of its possible states. It is therefore obvious that, despite the 
considerable differences in the formulations used, in either case one can speak about 
two possible ways of exhibiting one and the same fundamental property of quantum- 
mechanical objects. 

On the other hand, inequality (2) holds not only for density matrix, but as well 
for correlation functions of wavefields of various nature and it can be seminally applied, 
after appropriate terminology corrections, in other fields of physics. In fact, the basic 
works [1-3], where one-dimensional case of relation (2) was obtained, deal with partially 
coherent light beams. 

Given in [5] variant of the proof was based on the Carlson type integral inequalities 
for the Wigner function, that implies some disadvantages, in particular, because such 
approach does not determine the explicit form of the density matrix minimizing (2). 
This question is of interest in itself and, at the same time, is important for applications 
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(for example, at the coherence theory [8]). Hence in the present paper we shall consider 
an alternative way to derive the relation (2). The aim is not only to solve the problem 
stated, but also reinforces the inequality by taking into account the dependence of 
packing coefficient on the effective number of states C(s) =>• C(s,N e jj). 



2. Formulation of the problem and rigorous solution 

The method to be applied constitutes a modification of approach, which was developed 
in Ref. 9 to analyze an analogous one-dimensional problem. It is based on the use 
of standard Lagrange procedure of search for the minimum of the system uncertainty 
volume (i. e. the value AxAq) under a constraint of given N e g (5). As common for 
such treatments [1-3, 5, 9], one should seek the extremum of auxiliary functional F 

F[p(X,X')] = ^ 2 (Ax) 2 + l(A ? ) 2 , (8) 

which has the property that its minimum for k (k is a variate scale factor) is attained 
simultaneously with the minimum of uncertainty volume 

minF = 2 AxAq. 

The substitution of definitions (3), (4) into (8) shows that the value of F coincides in 
form with the mean value of energy of s- dimensional symmetrical harmonic oscillator 
with a Hamiltonian 

which is in the mixed state p(X, X'). It is therefore natural for further treatment to 
represent p(X, X') as a series expansion in the basis of this Hamiltonian eigenfunctions 

p(X,X) = ^a n , n ^ n (X)* n ,(X'), (9) 

n,n' 

where n = (ni, ...,n s ) is the "vector" index whose components range over nonnegative 
integer numbers, 



* n (X) = n^(a:i), 



i=i 



if> n (x) = J 2w J^p ex P {~k 2 x 2 /2) Hjkx), 

1 2 dj 1 2 

H n (x) = -^j=e, x ~r~^ e ~ X " ~ Hermite polynomials, the functions ^(X) are real and 
orthonormal. In view of the foregoing, the substitution of (9) into (8) gives 

F = ^a n , n f-||n|| + lY ||n|| = ^ni. (10) 
^ s / i=i 
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Now the additional condition, under which the effective number of states (5) is constant, 
can be written in the form 

A* = T?— = |a n ,n'| 2 = cons t (11) 

n,n' 

(in statistical optics the parameter /i is frequently referred to as global degree of 
coherence). One more constraint on expansion coefficients a njn / follows from the 
requirement for the density matrix (6) normalization 

Ev = 1 - ( 12 ) 

n 

Without going into details of elementary, although cumbersome intermediate 
calculations, we now turn to the analysis of the solution obtained. The first and most 
important consequence of minimization of F (10) in coefficients a n n is that in the state 
of least uncertainty p min (X., X') all off-diagonal elements of the matrix [a n ,n] are equal 
to zero 

^n,n' "n,n^n,n'' 

It means, that in this case the density operator's eigenfunctions [6] (in optics — the 
decomposition modes) coincide in form with the eigenfunctions of energy operator of 
multidimensional harmonic oscillator. Accordingly, the diagonal elements of matrix 
[a n ,n'] are the eigenvalues of density operator 

a n ,„^n(X) = J «'p mm (X,X / )^„(X / ). 

They define the probabilities to find the system in the pure state \l/ n (X) and satisfy the 
conditions Ima nn = (as a consequence of [a nj „/] is Hermitian) and < a nn < 1. 

It should be noted here that the use of characteristic N e g as a measure of number 
of possible pure states of the system is directly related to the operation of transition to 
the basis of density operator's eigenfunctions in the representation (9). In the theory of 
stochastic processes [10] this procedure is referred to as Karhunen-Loeve expansion of 
"correlation function" p(X, X') and gives, as known, the most compact and the most 
rapidly converging form of this series. Therefore, in this paper (as well as in [5]) under 
the term "state" is meant, as a rule, the eigenstate of density operator. 

The second conclusion, directly following from the form of (10) - (12), is that 
the weights of states a nn can only depend on the vector index norm ||n||, and this 
dependence is linear and decreasing when ||n|| increases. Then it is obvious that at any 
finite value of N e g the number of pure states Af present in the expansion p min with 
a probability other then zero is also finite. And in the index domain, the coefficients 
a n,n 7^ fill, layer by layer, the interior of s-dimensional equilateral pyramid with the 
total number of layers equal to C (0 < ||n|| < C). The most probable pure state 
(n = 0) corresponds to the pyramids vertex. The terms of expansion (9) pertaining 
to some particular layer are absolutely equivalent and make an equal contribution to 
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functional F. Hence the spectrum of the eigenvalues in this mixed state is degenerate, 
and degeneration multiplicity in each layer is determined by value of ||n|| and space 
dimensionality s 

The total number of significant terms in (9) is given by relation 

C— 1 

M(C) = £ g .(||n||) = = - S 9M- 

||n||=0 1 ' 

Finally, the sought state p m j„, which realizes the minimum of the generalized uncertainty 
relation, is described by series (9) with coefficients 



a n ,n = < M{C) 



i + (( £ -i).-IHI(. + i))^ (£) - JV * )( * + 2) 



,0 < || n || < C 



y/N eff s(C + s)(C-l) 

0, ||n|| > C 

(13) 

and the inequality analogous to (2) takes the form 

AxAq > i2C ; ' - " - W - ES = B CiN, s)) . (14) 

2(s+l) ( S + l)V7V^J(7+2) y 

And it should be kept in mind that quantity C, a certain positive integer, still remains 
a free parameter of the task and should be chosen proceeding from the condition of 
minimality of the right-hand side of inequality (14). By this means C will be a certain 
function of N e g. The above requirements for the expansion coefficients a n n (13) to be 
real and positive impose restrictions on allowable values of this parameter, namely, the 
quantity C should satisfy the following inequalities 

(C + s-l)\ s + 2 
eff > (£-2)!(s + l)!s + 2(£-l)" { } 

There exist certain ranges of values of N e g and s (in particular, s ^> 1), in which 
inequalities (15), (16) determine C uniquely. But where it is not fulfilled and several 
integers fall within the interval given by formulas (15) and (16), the sought value of C 
turns to be largest of them (i. e. integer proximate to the upper boundary of (16)). 

Inequality (14) is more hard than the previous inequality (2). On the plane of 
parameters with coordinates AxAg and N e ff it rigorously defines the region of physically 
realizable states. It is therefore apparent and can easily be proved that one should 
arrive at the same distributions of expansion coefficients (13) and the same region (14) 
by solving the inverse problem of seeking a mixed state with a largest possible value of 
A^g^ at a given measure Av p of phase volume occupied by this state. 
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3. Approximate form of uncertainty relation 

By virtue of the fact that with increasing N e g the quantity £ goes through a set 
of discrete values, the boundary of physical region B(N e jj) in (14) is not a perfectly 
smooth curve and is not described by analytical expression. Even though the numerical 
calculation of the right-hand side of inequality (14) presents no problem, it would be 
desirable to have its approximate analytical form in order to analyse the obtained 
relation and compare it with (2). 

To this end we have investigated the behaviour of the function B(N e ff, £) at a given 
value of N e ff, formally regarding £ as an independent continuous variable. Analytical 
and numerical calculations show that in the vicinity of the true value of £ the dependence 
of B on £ is extremely weak. As well, in this region slightly different neighbouring local 
minimum and maximum for B(C) occur. This gives grounds to use the continuous 
parameter £ instead of discrete number £ when approximately describing the boundary 
of inequality (14). There are several ways of choosing the specific condition determining 
the value of £. 

The following one seems to be fairly simple and logically substantiated. In 
the expression for expansion coefficients (13), one can formally require "continuous" 
transition of its first part to the second one|| which, after a little algebraic manipulations, 
yields a transcendental equation for determining £ in the form 

(s + 2)T(£ + s + l) 

a? p = =>■ M eff = = — ; 17 

£ ' £ ° (s + 2£)(s + l)!r(£)' V ; 

(because of £ being continuous the factorials entering into the formula for M{£) have 
been replaced here by Euler gamma-functions). Accordingly, with such a way of defining 
£, the approximate expression for inequality (14) can be written in the parametric form 

. . s + 2£(N eff ) , , 

H ~ 2(s + 2) y J 

It should be noted that the approximate boundary of physical domain, given by relations 
(17) and (18), exactly coincides with the result of the procedure proposed in [9], 
proceeding from the definition of £ as the upper bound of inequality (16). 

Another way of finding a smooth form of the boundary, B(N e ff), can be chosen by 
substituting £ by the value of £' corresponding to min^ B(N e ff, £). The thus obtained 
curve has the advantage that it is wholly located in the physical domain and is tangent 
to the strict boundary at individual points. Unfortunately, at large s this method does 
not lead to analytical formulas. 

Finally, according to the above mentioned property of equivalence of the direct and 
the inverse problems, the domain of physically realizable states can be represented as 
inequality to the largest possible value of N e ff at given AxAq. On the basic of (17), 

§ Naturally, in this case, the variable ||n|| in (13) should also be regarded as continuous. 
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(18) it gives the analytical expression 

Nff< 1 r(( 3 + 2)AsAg + a /2 + l) 
eff ~ 2AxAq (s + l)!r ((s + 2) AxAq - s/2) ' 1 1 

The introduction of the concept of "packing coefficient" C(s,N e ff) (whose explicit 
form is not described by simple analytical formula) permits to rearrange relations (14) 
and (18) into the form similar to (2). The point of such transformation is that thus 
the main dependence of generalized uncertainty relation on problems parameters is 
emphasized whereas C(s,N e ff) plays role of correction factor with comparatively weak 
dependence on N e g. 

The method of obtaining the approximate expression for the physical domain 
boundary supposes itself that the approximation (18) should tend to the exact formula 
(14) with increasing C (and, thus N e jj as well). Really, exactly in the limit of C ^> 1 
the substitution of a discrete value of £ by a continuous value of C slightly effects the 
weights of individual states in (13). However, even at small N e g > 1 the approximate 
relations (18) and (19) turn to be in a very good accord with the rigorous inequality 
(14). The reason for this is the above mentioned local behavior of B(C). 

The said is illustrated by plots of exact and approximate dependencies of C(s, N e g) 
given in the Fig. 1 for some values of s. From this figure it is seen that insignificant 
discrepancies between (14) and (18) only take place near N e jj pa 1 and they are the more 
pronounced, the greater the problem dimension s. Accordingly, as s is increased, thus 
approximate relation approaches the rigorous one evenly closely at larger values of N e g. 

Comparing inequalities (2) and (18), it can be easily shown that (2) is an asymptotic 
form of (18) (and, consequently, (14)) at N e g — > oo, but at s > 1 a good approximation 
of (2) to (14) is only attained at large N e g. At the same time, the refined form of the 
uncertainty relation (14) correctly describes all range of values of N e g and, in particular, 
the ultimate case of pure states 



C(s,N eff = 1) = 1, (AxAq) 



> 

N„„=l 



The coincidence of the asymptotic form of packing coefficient C(s,N e jj ^> 1) with 
formula (7) justifies the use in [5], when proving (2), of assumption that the Wigner 
function is nonnegative for the density matrix in the state with minimal uncertainty. 

4. Discussion 

At a qualitative level, the treatment of standard uncertainty relation as a requirement 
for phase space quantization in wave mechanics is rather frequently occurred in the 
literature on physics (see, for example [11]). Therefore, it can serve, in some way, as an 
argument in substantiating the postulate of statistical mechanics on the number of cells 
in phase space of the system: 

^ = W, (20) 
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where V P (E) is the phase volume of the system with given energy E and iV is the 
number of levels with energies not exceeding E (detailed definitions can be found in [4]). 
The results of [1-3, 5, 9] and the present paper open up the possibility of quantitative 
comparison of relations (2), (14) and (20). 

First of all their similarity lies in the fact that the quantities of the same physical 
nature enter into the right- and the left-hand sides of expressions (2) and (20) by identical 
manner. But the differences between (2) and (20) are far more essential and there is 
nothing strange about it, if we take into account the way in which the concepts of the 
phase volume and the number of states are introduced in both cases. 

The postulate given by (20) specifies the relation between energetic structure of 
the quantum system and the phase space which can be associated with this system in 
the quasi-classical approximation [4]. Such a relation is assumed to be universal and 
independent of a particular Hamiltonian of the system. Formula (20) describes system 
as a whole and has no connection to its particular physical state. 

On the contrary, inequalities (2), (14) refer just to the physical state of quantum 
system. Like Eq. (20), they are explicitly associated with neither the kind of 
Hamiltonian nor the structure of energy levels. The only limitation is that since the basic 
functions of expansion (9) form a complete set, the wave-packets, satisfying the condition 
of minimum uncertainty for mixed states, can be constructed only for Hamiltonians 
whose eigenfunctions also form a complete orthonormal basis. This holds, in particular, 
for the simplest system of noninteracting particles in a free space. In general, however, 
the eigenstates of density operator are not the states with constant energy and, thus, 
the state with minimum uncertainty will not be stationary. The exception is the case 
of multidimensional harmonic oscillator. 

It is clear from the said that the quantities Av p and V p , respectively, in Eq. (2), 
(14) and (20) not only do not coincide numerically with one anther (even for oscillator), 
but also have a different operational meaning. A similar statement in general case is 
true for the quantities N e g and N as well. So there is nothing surprising in the fact 
that relation (2) has the form of inequality, while (20) — equality. Possibly, a better 
correlation between (2) and (20) can be attained by using other, alternative formulations 
of the uncertainty principle [9]. 

It would be worth to draw attention to one more difference between (2) and (20). 
In the postulate given by (20), the phase space volume of any quantum cell is constant 
and varies with change of dimensionality in a strictly definite way — as (const) 5 . On the 
contrary, in the generalized uncertainty relation (14), the minimal specific phase volume 
(i. e. the volume per one effective pure state of the system) depends on the value of 
N e g. As seen from Fig. 1, the quantity C(s,N e jj) is less than unity and monotonically 
decreases with increasing N e g. It means that as the number of pure states involved 
in the statistical mixture for p(X, X') raises, there is a gain in the packing density of 
states. 

In the paper [5] it was also noted the effect of increasing packing with augment of 
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the configurational space dimension s, that is expressed by the inequality 

C(ks) < C(s) k , (k > 1). 

From (19) one can easily see that an analogous relation for the coefficients C(s, N e g) in 
the general case does not hold and this property of packing is asymptotic, i. e. in each 
particular case (s, k) it takes place beginning with sufficiently large N e g. The reason 
for this should be sought in the peculiarities of behavior of degeneration multiplicity 
<7 s (||n||) as a function of ||n|| and s. It is natural to suppose that the decrease in the 
coefficient C(s) with increasing s is due to the rapid growth of the value of g s , and for 
different dimensions it would be proper to compare the packing at equal values of the 
parameter £ rather than N e g. 

In this context, it is interesting to consider the case of the density matrix (9) with 
equal weights of all pure states a n n / = 5 n;n '/Af involved, that somewhat resembles the 
definition of iV in (20). For this situation all methods of determining N e jj gives N e g = Af 
and, although the minimum of uncertainty relation (14) for such a system is slightly 
exceeded, both the above properties of the packing coefficient are fulfilled. 

In conclusion, it is necessary to discuss also the question of what physical meaning 
is attributed to the arguments (coordinates) entering into the density matrix p(X,X'). 
These may be either the coordinates of one particle in the real three-dimensional 
space or the coordinates of several particles in the configurational space or both 
at once. The above solution for the density matrix in the state with minimal 
uncertainty is invariant under any permutation of its arguments (x, ^ Xj). For the 
real coordinates of one particle it means that the corresponding wave-packet possesses 
some rotational symmetry in three-dimensional space. But when arguments being 
interchanged correspond to several particles it is necessary to take into account the 
properties of their permutational symmetry. For indistinguishable identical particles 
the state with such a symmetry can only be realized in the case when these particles 
follow the Bose statistics. On the contrary, identical particles with Fermi statistics can 
not be described by a density matrix of the form of (9) with the coefficients given by 
(13). The state with minimal uncertainty for fermions should be sought from the very 
beginning in the class of wavefunctions antisymmetrized in permutations of arguments, 
which, naturally, should lead to a result drastically different from (14). 
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Figure caption 

Fig. 1. The packing coefficient C(s,N e jj) vs. the effective number of pure states 
in statistical mixture (9) for dimensions s = 1, 2, 3. Solid line — strict inequality 
(14), dotted line — approximation (18), dashed line — the asymptotic value C(s) for 

Neff » 1 (7). 



